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We study the time evolution of ultra-cold atoms in an accelerated optical lattice. For a Bose- 
Einstein condensate with a narrow quasi-momentum distribution in a shallow optical lattice the 
decay of the survival probability in the ground band has a step-like structure. In this regime we 
establish a connection between the wave function renormalization parameter Z introduced in [Phys. 
Rev. Lett. 86, 2699 (2001)] to characterize non-exponential decay and the phenomenon of resonantly 
enhanced tunneling, where the decay rate is peaked for particular values of the lattice depth and 
the accelerating force. 



PACS numbers: 03.65.Xp, 03.75. Lm 

I. INTRODUCTION 

Resonantly enhanced tunneling (RET) is a quantum 
effect in which the probability for the tunneling of a par- 
ticle between two potential wells is increased when the 
quantized energies of the initial and final states of the 
process coincide. In spite of the fundamental nature of 
this effect p] and its practical interest [2], it has been 
difficult to observe it experimentally in solid state struc- 
tures. Since the 1970s, much progress has been made in 
constructing solid state systems such as superlattices j3J- 
[5] and quantum wells [6] which enable the controlled ob- 
servation of RET [7] . 

In recent years, ultra-cold atoms in optical lattices [HI 
9 , arising from the interference pattern of two or more 
intersecting laser beams, have been increasingly used to 
simulate solid state systems [9"rfTT]. Optical lattices are 
easy to realize in the laboratory, and the parameters 
of the resulting one-, two- or three-dimensional periodic 
potentials (the lattice spacing and the potential depth) 
can be perfectly controlled both statically and dynami- 
cally. In the papers [HI [12] , a Bose-Einstein condensate 
(BEC) in accelerated optical lattice potentials was used 
to study the phenomenon of RET. In a tilted periodic 
potential, atoms can escape by tunneling to the contin- 
uum via higher-lying levels. Within the RET process the 
tunneling of atoms out of a tilted lattice is resonantly en- 
hanced when the energy difference between lattice wells 
matches the distance between the energy levels in the 
wells. 

The atomic temporal evolution is described by the 
survival probability, starting from an initial state pre- 
pared in the ground band of the lattice. At long interac- 
tion times, after several tunneling processes, the survival 



probability is characterized by an exponential decay rate 
with a constant tunneling probability for each Bloch pe- 
riod [14] - Such a decay was examined in different theo- 
retical analyses [S[ [H] and measured in experimental 
investigations with ultra-cold atoms [12 [T31 EH US] . In 
this study we scrutinize the time behavior of the tunnel- 
ing probability and use its remarkable features at short 
and intermediate times in order to extract information 
about wave-function renormalization effects. 

The key quantity in this context is the probability that 
the system investigated "survives" in a given state (or a 
set of states, such as a band of a lattice). In this article 
we shall deal with survival probabilities whose behavior 
is complex and difficult to analyze. See for example the 
experimental results of ref. [TS] and the Figs. [2] and [5] in 
the following, which display the survival probability of a 
cloud of ultra-cold atoms in the ground band of an ac- 
celerated optical lattice. Clearly, one can properly speak 
of the "decay" associated with an unstable system (the 
atoms tend to leak out of the accelerated lattice) , but the 
time evolution can display oscillations or even plateaus. 
(As we shall see, the latter are easily understood in terms 
of the initial atomic state.) 

General theoretical considerations show that the (adi- 
abatic) survival probability of an unstable system can 
often be written as 

P(t) = Z exp (— 7f) + additional contributions, (1) 

where 7 is the decay rate, which can be computed by the 
Fermi golden rule, and the parameter Z, representing the 
extrapolation of the asymptotic decay law back to t = 0, 
is related to wave-function renormalization. Law ([!]) is 
valid both in quantum mechanics [17j [18] and quantum 
field theory [THl [20] , and Z can be smaller or larger than 
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unity [3T|. Typically, the additional contributions in ([!]) 
dominate both at short and long times, where the expo- 
nential decay law is superseded by a quadratic [2"2"H2"i] 
and a power law [55], respectively. They are therefore 
crucial in order to cancel the exponential in these time 
domains. However, they can play a key role in a much 
more general context, such as the RET phenomenon to 
be investigated in this article. 

The pioneering experiments performed in Texas, with 
Landau-Zener transitions in cold atoms, checked the ex- 
istence of the short-time quadratic behavior |26| and the 
transition [27 from the quantum Zeno effect [53] to the 
anti- or inverse-Zeno effect [58H5U] . through a sequence 
of properly tailored quantum measurements. 

With the arrival of Bose-Einstein condensates the ex- 
perimental resolution has advanced even further as com- 
pared to cold atoms. While cold atoms can have a mo- 
mentum distribution on the order of a Brillouin zone or 
more, a very narrow distribution (much smaller than a 
Brillouin zone) is achievable with BECs. Even the step- 
like structure of the survival probability occuring for shal- 
low lattice depth can be resolved with great precision 
[15] [T6] . It is in this regime of shallow lattices and short 
jump times [31] where the yet unobserved link of RET 
and the initial deviation from exponential decay is most 
striking. This work is devoted to the study of these ef- 
fects. The choice of a different initial atomic state, with a 
well defined momentum, will enable us to observe a more 
complicated temporal structure. We shall therefore scru- 
tinize the time evolution in order to unveil an exponen- 
tial regime and introduce the Z parameter in our RET 
framework. 

The paper is organized as follows. We briefly sum up 
previous results on RET and the quantum Zeno effect in 
Section [TTj We then analyze the dynamics in the tilted 
lattice in Section III and show in Section IV the main 



part of this article, how the two phenomena arise as in- 
terference effects. Section [V] reports experimental results 
for the wave-function renormalization parameter Z in the 
case of a Bose-Einstein condensate in an accelerated opti- 
cal lattice, and also a comparison with the experimental 
configuration by Wilkinson et al. [55]. Section VI con- 
cludes our work. 



II. LANDAU-ZENER AND RESONANTLY 
ENHANCED TUNNELING 

A Landau-Zener (LZ) transition takes place in a sys- 
tem with a time-dependent Hamiltonian, in which the 
spectrum, as a function of a control parameter (here time 
t), is characterized by the presence of an avoided cross- 
ing [35-36] . A LZ transition is described by the following 
two-level Hamiltonian 



written in a suitable basis, known as diabatic basis. The 
expectation values of Eq. ([2| on the two states of the basis 
depend linearly on time and cross at t = 0. On the other 
hand, the coupling SE/2 between the states is constant. 
The diagonalization of Eq. ([2| yields the eigenvalues 



E ± = ±\ (at) 2 + 



SE 



(3) 



The eigenbasis of H^z(t) is called the adiabatic basis. At 
t — > — oo the adiabatic energy levels of Eq. {3} are in- 
finitely separated, and no transition between them oc- 
curs. The distance between the levels decreases towards 
the avoided crossing at t = 0, and then increases again 
until, at t — > oo, the separation becomes again infinite. If 
the system is prepared at t — > — oo in one of the adiabatic 
eigenstates, the probability that the system undergoes a 
transition at t — > oo towards the other adiabatic eigen- 
state reads EH 



P L z = exp - 



n(SEf 
Aha 



(4) 



( at 6E/2 \ 
HlzW- \ §e/2 _ at I 



(2) 



A particle in a shallow periodic potential, subjected to 
an external force, is an example of a physical system in 
which a LZ process can be observed. In this case, the di- 
abatic basis is represented by the momentum eigenstates. 
As schematized in Fig. [T] if the system is initially pre- 
pared in the lowest band, with a very peaked momentum 
distribution around p = 0, it will evolve towards the edge 
of the first Brillouin zone, where the distance between the 
first and the second band is minimal and transitions are 
more likely to occur, and then evolves back to the bottom 
of the first band. The transition probability towards the 
second band in this process can be approximated by Plz, 
but discrepancies can arise due to the differences between 
the idealized case, leading to the LZ formula Q , and the 
real physical situation. Indeed, the periodicity of the lat- 
tice implies that the aforementioned process occurs in a 
finite time, and that in the initial and final states the 
adiabatic levels are not infinitely separated. The correc- 
tions to the LZ transition probability due to the finite 
duration of the process are discussed in [16l [37] . 

Other corrections to Eq. Q should be considered if the 
lattice is not shallow. In this case, couplings to higher 
momentum states play an important role and a two-level 
description is not a good approximation anymore. 

Moreover, there is another kind of deviation from LZ, 
which will be the main object of our analysis. Since Eq. 
Q is obtained under the hypothesis that only one of the 
two adiabatic eigenstates is initially populated, it is not 
valid anymore if both states are populated. These devia- 
tions can be relevant even if one of the initial populations 
is very close to zero, since their order is square root of the 
smaller population, as will be discussed in the following. 
In a periodic potential, tilted by an external force F, the 
probability that a wave packet initially prepared in the 
first band jumps to the second band corresponds to the 
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FIG. 1: (color online) Energy diagram for a particle in a pe- 
riodic potential vs either the time, in units of Bloch time Tb 
defined in (14 1, and the atomic momentum p a t/2p r ec where 
free = hn/d,L units. Under the application of an external 
force, the quasi-momentum increases with time and at the 
avoided crossing between two bands, at the edge of the Bril- 
louin zone, where p a t = Prcc, the condensate tunnels with a 
probability amplitude ptj and survives with an amplitude Sij, 



defined in Eq. ( 22 1 below. Between two avoided crossings of 
the lowest two bands a relative phase (j>, defined in Eq. (17 1, 
is acquired, which is graphically displayed in the figure as a 
grey area. The final survival probability at a given time is 
the sum over all possible routes, just like a path integral in 
momentum space. 



LZ prediction Q only if the second band is empty. A 
small population in the second band gives rise to oscilla- 
tions around Plz- 

Finally, the transition probability is enhanced by a 
large factor with respect to the LZ prediction if the en- 
ergy difference Fc^Ai between two potential wells (d,L 
being the lattice spacing and d^Aj the distance between 
the wells) matches the average band gap of the non- 
tilted system (RET). One expects that in a RET process 
from the first to the second band, the asymptotic regime 
will only be reached after a transient period. Indeed, 
while the first transition occurs when the second band 
is strictly empty (and thus the tunneling event closely 
follows the LZ prediction), further RET transitions will 
occur periodically in time and, starting from the second 
tunneling process, interference effects due to the finite 
population amplitude in the second band will start to 
play an important role, modifying the time evolution in 
an important way. 

The analysis of the following two sections will endeavor 
to take all these effects into account. We shall build 



up an effective model, whose validity will be tested for 
rather diverse ranges of the parameters, and compared 
to experimental results finally in section W\ 



III. DYNAMICS OF INTERBAND TUNNELING 

In our analysis we are interested in describing the RET 
process from the first to the second band of a Bose- 
Einstcin condensate loaded into an optical lattice. It will 
be assumed that almost all the particles of the system 
are in the condensate, so that the system is described by 
a single-particle wave function ip(x,t) [38 . Moreover, we 
will consider the condensate dilute enough so that the 
interaction between particles can be neglected. This im- 
plies that the wave function of the system obeys a linear 
Schrodinger equation. Nonlinear effects have been stud- 
ied in the RET regime in previous works [12 US l3"9TEH] . 

The experimental condition is that of an accelerating 
one-dimensional optical lattice, with constant accelera- 
tion a. In the rest frame of the lattice, a particle of mass 
m sitting in the lattice is subjected to an external force 
F = ma, and thus the time-independent Hamiltonian of 
the system in this frame of reference reads 



H = — 



h 2 d 2 v (2-kx\ TT _ 

7, — jr^ + 7T cos ~3 — ) — Fx = H — Fx, 
2m ox z 2 \ ul J 



(5) 



where V is the lattice depth and the lattice period d^ 
is half wavelength of the counterpropagating lasers. Hq 
represents the "unperturbed" Hamiltonian, whose eigen- 
states are the Bloch functions 



<t>a,k{x) = e ikx x a ,k(x), 

Xa,k{ x + d L ) = Xa,k{x), 

H Q (f> atk (x) = E a (k)</> at k(x), 



(6) 
(7) 
(8) 



with a — 1,2,3,... the band index and k the quasi- 
momentum, ranging in the first Brillouin zone B = {k\ — 
tt/c^l < k < it jdh). The dynamics of the system depends 
on two dimensionless parameters |16j . related to lattice 
depth and external force: 



V 



V 



E„ 



,F 



Fd L 

Ey r 



with E Te 



2m 



7T 



(9) 



where m is the mass of the atoms. Applying the unitary 
transformation 



i/j(x,t) = exp(—iFxt/h)ip(x,t) 



(10) 



restores the translational invariance of the Hamiltonian, 
at the expense of an explicit time dependence: 



H(t) 



1 

2m 



-ih 



d_ 
Ox 



Ft 



V 



{2irx\ 

— cos — - — 
2 \d L ) 



(11) 



Rewriting this Hamiltonian in the momentum basis as in 
[16] establishes the relation to the Landau-Zener Hamil- 
tonian introduced in Eq. To calculate the time evo- 
lution of any momentum eigenstate, we only need the 
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Hamiltonian Hk acting on the subspace with a given 
quasi-momentum ha, as there is no transition between 
states with different k 



( 



2m 







'(*-£) a 
mV/2 



mV/2 

(hk) 2 

mV/2 



mV/2 



V o 



■■■/ 

_ (12) 

where k = k(t) = k n + Ft/h. This Hamiltonian (12) leads 



to a very accurate numerical solution of the Schrodinger 
equation. For small V on the scale of E rec its dynamics 
is well described by successive Landau-Zener transitions, 
occuring whenever two diagonal terms in H^ become 
degenerate. We will use this approximation to obtain 
analytical results. In Fig. [I] the relevant transitions are 
depicted graphically. 

We first examine an adiabatic approximation of the 
dynamics generated by the Hamiltonian d5j), yielding no 
transition between bands (single-band approximation), 
which will highlight the time periodicity of the system 
and the phase differences between bands. We shall 
then introduce an effective coupling between the low- 
lying bands that enables one to obtain transition rates. 
The adiabatic approximation is consistent if Fd^ < V, 
namely if F Q < V in Eq. 

The initial state will be assumed to be highly peaked 
around a single quasi-momentum value kg, that is, the 
width of the initial quasi-momentum distribution will be 
taken to be much smaller than the width 2ir/dL of the 
first Brillouin zone B. In this situation, it can be proved 
[361 137] that in the adiabatic single-band approximation 
the average quasi-momentum evolves semiclassically, so 
that at time t 



k(t) = k + 



Ft 



(13) 



with negligible spread in the quasi-momentum distribu- 
tion occurring during the evolution. This yields Bloch 
oscillations in a tilted lattice with a Bloch period 



2irh 
FaW, 



h 2tt 



(14) 



The initial state analyzed here has a well defined initial 
momentum (in 23), but can be distributed among differ- 
ent bands. At the end of each Bloch period, the ampli- 
tude in band a acquires the following phase with respect 
to the amplitude in band /3 



A<p al3 = -~(E a (k) - E p {h))Ta = 



2tt 



,(*) 



(15) 

where (...) denotes the average over B and e 7 (fc) is the 
energy of the state with quasi-momentum k in band 7 in 
units -E rec - 

We now analyze inter-band transitions through an ef- 
fective model. We focus on the experimental parameters 
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FIG. 2: (color online) Adiabatic survival probability in the 
lowest band P(t) vs. time, obtained by numerically solving 
the Schrodinger equation for the atomic evolution under the 
Hamiltonian of Eq. (111. The initial state has initial quasi- 
momentum 0.2p rcc in the Brillouin zone and a negligible quasi- 
momentum width. Upper panel: Vb = 1.5,0 ~ 47r; lower 
panel: Vb = 4.5, <f> ~ 2-k. The solid (blue) line represents the 
unperturbed time evolution, the dashed (green) line the time 
evolution with the phase change discussed in Sec. |V B| after 
each Bloch period. The red crosses are an extrapolation of 
the first time step to the following periods. Plateaus manifest 
for Vb = 1.5 (upper panel). The validity of this picture is 
discussed in the text and in Appendix [X] 



of the Pisa setup [HI [16] and model transitions from the 
first to the second band. In the parameter regime of shal- 
low lattices there is numerical and experimental evidence 
of a step-like structure of the adiabatic survival proba- 
bility P(t) [H] in the first band. If the initial state is 
peaked around k = and lies in the first band, the sur- 
vival probability is characterized by steep drops around 
times t — Tb{u + 1/2) with n integer, and flat plateaus 
between these times [40]. This view is corroborated by 
numerical simulations (Fig. [2j upper panel) and experi- 
mental observations 1151. This time evolution is due to 
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the fact that the coupling between the first and the sec- 
ond band is maximal at the edge of the first Brillouin 
zone, for k = dzn/d^, and thus significant transitions 
occur there, with periodicity Tb. Figure [2] shows that 
plateaus are clearly present for Vq — 1.5 (shallow lattice, 
upper panel), but start to wash out for Vo — 4.5 (lower 
panel) . The range of validity of the plateau picture is fur- 
ther discussed in Appendix [X] and is approximately valid 
for Vq < 4.5. In the following analysis we shall focus on 
this regime. 

The approximated dynamics takes into account exper- 
imental and numerical evidence and is valid for small 
values of Vo and Tb, when the transition times can be 
considered much smaller than Tb. We assume that the 
evolution inside the first band is adiabatic for all k, except 
for k ~ n/di,, when a transition towards the state with 
the same quasi-momentum in the second band becomes 
possible. This transition will be effectively described by 
the evolution operator of the form 



U 



Sl2 -Pl2 
Pl2 S12 



(16) 



with pi2 = yl — s 2 2 . The operator U acts on the two- 
dimensional space spanned by {|1) , |2)}, where |1) rep- 
resents the state with k = tt/cIl in the first band and 
1 2) the state with same quasi-momentum in the second 
band. 

The transition from the second to the third band can 
be schematized as the loss of a fraction 1 — s| 3 in the 
population of the second band towards a continuum, oc- 
curring at the crossing around k — 0. This assumption 
is justified for small values of Vb (see discussion above), 
such that a particle in the third (or higher) band can be 
considered free. 

During each Bloch cycle separating two successive 
transitions, the relative phase between the second and 
the first band amplitudes increases by Jl5|), which reads 



(17) 



cj>(V ,F ) = ^(AE(Vo)), 



where (AE) is the energy difference (in units T roc ) be- 
tween the second and the first band, averaged over B. 
This quantity can be exactly computed by using Math- 
ieu characteristics a(n,q), which are the eigenvalues of 
the Mathieu equation [35] 



d y 
dx 2 



+ [a - 2qcos(2x)] y = 0, 



(18) 



corresponding to the Floquet solutions y(x) — 
e-xjp(inx)u(x) . For small Vq, a good estimate is given 
by a Landau-Zener gap integration 



(AE) 



64 



V 2 



-arcsmh-. 



(19) 



For larger Vo, a tight-binding, or harmonic oscillator, ap- 
proximation yields 
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FIG. 3: (color online) Average band gap (AE) vs Vo, both in 
E ICC units. Com pari son between numerical results, analytical 
results from Eq. ( 19 1, and the harmonic oscillator approxima- 
tion (201. In the small V regime the band integration yields 



a good approximation, while for larger V, where the coupling 
becomes continuous, the harmonic oscillator approximation is 
more effective. 



The exact result and the two aforementioned approxima- 
tions are compared in Fig. [3j 

The effects of the dynamics in a time Tb from one 
transition to the next one can thus be modelled in the 
basis {|1) , 1 2)} by an effective non- unitary operator 



W 



1 







s 23 e l 



(21) 



By making use of this simplified model, we describe 
the time evolution in the following way. At t = the 
condensate is in the first band, with quasi-momentum 
close to k = 0. As the lattice is accelerated, the quasi- 
momentum increases until it reaches n/di, at t = Tb/2, 
where the operator U comes into play and transfers part 
of the population to the second band. The evolution from 
Tb/2 to 3Xb/2 is summarized by the application of W. 
Then, the second transition occurs, and part of the pop- 
ulation in the second band (decreased by losses towards 
the third band) can tunnel back to the first band due to 
the action of U and gives rise to interference effects. The 
same steps occur in the subsequent transitions. 

On a time span Tb, the dynamics of the system is 
therefore determined by the successive action of the non- 
unitary operator 



i<t> 



(AE) 



1. 



(20) 



U = UW = I - Jrz yr ^ 6 I (22) 
in the basis {|1) , |2)}. The order of the two operations is 
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not relevant, since W acts trivially on the "initial state" 
|1) before the first transition. 

Besides the phase </>, the operator U depends on two 
other independent parameters, namely the survival am- 
plitudes S12 and S23. S12 represent the survival amplitude 
in the first band after the first transition, which is in fact 
comparable to a LZ process since the second band is ini- 
tially empty. The survival probability S23 is related to a 
LZ tunneling from the second to the third band, if we as- 
sume the third band to be empty before each transition 
process. A graphical representation of the parameters 
appearing in Eq. (22) is given in Fig. [lj 



Using the LZ critical acceleration for the first and sec- 
ond band gap [32j [33j 06l 07] , analytical expressions for 
S12 and S23 as functions of the microscopic parameters 
can be obtained. At lowest order in Vb, the survival am- 
plitudes read 



S12(Vq,F ) 



S23(V ,F ) 



(i j) 

where P^z i s the Landau-Zener transition probability 
Q from band i to band j. 

The evolution on a timescale Tr, determined by a se- 
quence of U operations, will be analyzed in detail in the 
following section. 



l-P^\v ,F ) 



= Wl - exp - 



32F 



P& 3 \Vo,F ) 



1 - exp - 



32- 16 2 (2F ) 



(23) 



(24) 



IV. TRANSIENT AND ASYMPTOTIC 
BEHAVIOR 



1$. 



cie^i) +c 2 e 2 l |V'2) 



(27) 



Due to the dissipative term in W, the two eigenvalues 
are smaller than unity, and one of them, say e\, is larger 
in modulus than the other one. Thus, for n sufficiently 
large, the evolution reaches an asymptotic regime, in 
which the state after the n-th transition is determined 
only by the state after the previous one, with a transi- 
tion rate depending on the largest eigenvalue. Since the 
survival probability in the first band can be defined as 
P n = j (1|$„) | 2 , in the asymptotic regime one gets 



P n ~ |ei| 2 P„_i. 



(28) 



By defining an asymptotic transition rate 

7 = -log(| ei | 2 ), (29) 

it is possible to introduce a function P z (i) that coincides 
with the value of the survival probability at the center of 
the plateaus, at times t = nT^: 



P z (t) = Ze X p{-jt). 



(30) 



Compare with Eq. Q. The parameter Z in Eq. (30) 
is in general different from unity, due to the transient 
regime at the beginning of the evolution. It represents 
the extrapolation of the asymptotic exponential proba- 
bility back at t = 0. 

We now derive an analytical expression for Z. In the 
asymptotic regime, the system evolution described by 
Eq. (27) corresponds to an evolution operator applied 
to an initial unnormalized vector |^o) = c-i ^i)'- 



We now specialize the model outlined in Section III 
to the Pisa experimental setup [15] [16] . The state of the 
system before the first transition is |1). Immediately after 
the n-th transition, occurring at time t = T&{n + 1/2), 
the state of the system is 



|*n) =C/ n |l) 



(25) 



The matrix U in Eq. (22) can be diagonalized, yielding 



eigenvalues (ei,e2). By expanding the initial state as 

|l) = c 1 |^ 1 )+c 2 ^ 2 >, (26) 

where IV'1,2) are the normalized non-orthogonal eigenvec- 
tors of U, the state of the system at time T^(n + 1/2) 



|$„) ~ dej = U n ( Cl |Vi)) = U n |* ) 



(31) 



The Z parameter, representing the extrapolation of the 
asymptotic behavior back to t = 0, can be defined as the 
square modulus of the projection of the fictitious initial 
vector 1^0)1 onto the actual initial state |1) 



(l|*o)| 2 



(l|^i)| 2 



(32) 



which corresponds to an extrapolated "survival proba- 
bility" in the subspace spanned by |1), evaluated at the 
initial time. Z can be analytically computed as a function 
of the independent parameters of the model, by explicitly 
diagonalizing U. One obtains 
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with 



Z(S\2, S23, 



^(si2,s 23 , 
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^(1 - s 23 cos0) - 


. / K.(s 12 ,S 2 3,4 ) ) 


+ s 2 3 sin 2 (j> 


2— Sj 2 (l + S 2 3 COS0) 


1 ^12 




h V 8 


\/2/C(si2,S23,0) 
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K(si2.S23. 



'C(s 12 ,s 23 , 



[2 - s\ 2 (l + s 23 cose 



S 12 (1 + 2s 23 COS <f) + S 2 3 cos(2</>)) — 4s 2 3 cos ^ 

sf 2 (1 + 2S23 COS ^ + s 23) — 8s 23 S 2 2 (COS (j) + 2s 23 + s 23 cos ( 



16s 2 3 . 



(33) 



(34) 



In order to gain a qualitative understanding of the de- 
pendence of Z (and 7) on the phase difference 4> acquired 
during a Bloch cycle, let us compare the first and second 
transitions. Let P = 1 be the initial value of the survival 
probability in the first band. After the first transition, 
the survival probability becomes 



P a = s( 2 P = c 



'Pa 



(35) 



At the second transition, the discrepancy with the LZ 
prediction becomes manifest. Since, in the parameter 
regime of small Vq we are considering, the ratio S23/S12 
is very small [see Eqs. ( 23 )-( 24 1] , we can apply a first- 



order approximation, yielding 



P2 



(S 2 2 - 2s 2 3Pl2 C0S ' 



'Pi- 



rn 



Thus, if the phase is (f> = 2nj, with j E Z, the second 
transition is enhanced with respect to the first one. In 
this local maximum in the transition rate as a 

function of Fq is expected. On the contrary, if <\> = (2j + 
l)7r, the second transition is less pronounced than the 
first one. 

A backwards extrapolation of the second step gives a 
rough estimate of the Z parameter, which we call Z\: 



Z x 



1 + 2s 23 



P12 
S12 



cos < 



(37) 



Even if Eq. ( 37 1 represents a rather crude approximation, 
it brings to light the correspondence between resonances 
in the asymptotic transition rates and resonances in the 
Z parameter. Quantities like (37) are very useful in an 
experimental context, where only the first few steps in the 
Bloch cycles are accessible. If the survival amplitude can 
be measured up to the ./V-th transition, the Z parameter 
can be approximated by 



Z 



N 



At the same time, 



7 — In- 



(38) 



(39) 



The convergence to the real value of Z is typically very 
fast, and the first few cycles are already sufficient to ob- 
tain an excellent approximation. 



The estimates of Eqs. (17|-(23|, together with Eq. (33) 



enable one to obtain an analytical expression Z(Vn, Fq), 
yielding the value of Z as a function of the microscopic 
parameters. Figure [4] shows a comparison of the numer- 
ical calculation and the estimates for 7 and Z with our 
analytical model. It is clear that the model yields a bet- 
ter approximation for smaller Vq. For Vq ^ 4.5 the peaks 
of Z are overestimated and the picture of successive tun- 
neling events with an intermediate phase accumulation 
becomes less valid. In the regime of small Vq, the analyt- 
ical model is very efficient, as long as Fq is not too large 
and the LZ tunneling rates do not have to be adjusted 
due to the finite initial time of the evolution [37] . 



V. EXPERIMENTAL CONFIGURATIONS 

This section contains a discussion of the experiments 
performed up to now and suggestions for future mea- 
surements aimed at controlling the decay by a manipula- 
tion of the phase of the temporally evolved atomic wave 
packet. The relations of Sec. |IV| can be tested experi- 
mentally as follows. 



A. Measurement of P(t) 

An experimental check of the theory at the basis of 
the wave-function renormalization Z is obtained by mea- 
suring the survival probability P(t) for a time up to five 
Bloch periods for different parameter values, as in Fig. 
[2] and then introducing a fit with the exponential law of 
Eq. ( 30 1 for the survival probability at times t — nT^ ■ 



The Z and 7 parameters are determined by such a fit. 
The results of this approach are discussed in the following 
for the case of a narrow atomic momentum distribution, 
as in the RET experiments at Pisa with a Bosc-Einstein 
condensate [12l Q31 |T5l [16] , and for the case of a broad 
atomic momentum distribution, as for the experiment 
performed at Austin [57] . 
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1. Pisa RET experiment 

The time dependence of the adiabatic survival prob- 
ability was measured by freezing the tunneling process 
through projective quantum measurements on the states 
of the adiabatic Hamiltonian [15] . Experimental results 
of P{t) for different values of the lattice depth and the ap- 
plied force are shown in Fig.[5j The solid and dashed lines 
are a numerical simulation of our experimental protocol 
and an exponential decay fit for our system's parameters, 
respectively. The vertical intercept of the exponential de- 
cay at t = gives the value of and the exponential decay 
rate gives the value of 7. 

The resonant tunneling appears as a strong variation for 
the exponential decay rate of 7 as a function of </>, as mea- 
sured in the experiments |12[ 1 1 3 j . This variation matches 
the numerical predictions of Fig. |4| 

Measured values of the Z parameter vs the <f> param- 
eter are plotted in Fig. [6]ja). The error bars on the Z 
values are determined by the exponential fits, as in Fig. 
[5] Notice that Z values both larger and smaller than one 
are measured. The error of the phase <f> is linked to the 
experimental accuracy of the Vq and Fq parameters (Vq 
carries an error of around ten percent). The experimen- 
tal results are compared to theoretical predictions for the 
numerical solutions of the time-dependent adiabatic sur- 
vival probability. The peaks in the plot are determined 
by RET resonances. The simulation of Fig. [4] evidences 
that the dependence of Z on <j> matches the dependence 
of 7. The position of the largest peak corresponds to the 
main resonance [T2l [T3] Ai — 1, and the positions of the 
smaller peaks are in agreement with those of higher order 
resonances. The agreement between the theoretical and 
experimental determinations of Z is very good, taking in 
account the difficulties of a precise determination of the 
lattice depth Vq. It should be noticed a posteriori that 
the experimental results are more easily produced in the 
case Z < 1. 



2. Austin experiment on non-exponential decay 




1.5 2 2.5 
< A E > /F 



FIG. 4: (color online) Decay rate 7 and wave-function renor- 
malization Z vs {AE)/Fo = cj>/2n. Comparison among an- 
alytical results, obtained by exact diagonalization of the re- 
duced evolution operator U in Eq. (22 1 [(red) dashed lines], 
numerical simulations based on Eq. (Ill (solid lines) and 



The very broad atomic distribution of the experiment 
perfomed by Raizen's group in Texas [26l EI], occupy- 
ing several Brillouin zones, leads to a different temporal 
evolution of the survival probability. In particular, the 
deeper lattice potentials used in these works imply a dif- 
ferent behavior of the Z function. The survival proba- 
bility was numerically evaluated on the basis of the the- 
oretical treatment reported in Niu and Raizen |48) and 
Wilkinson et al. For the case of Rb atoms and pa- 
rameters very close to those experimentally investigated 
in Pisa, Fig. |6jb) reports the Z function versus the pa- 
rameter at a fixed value of the lattice depth. It may 
be noticed that the values of \Z — 1| are smaller than 
those measured in the case of a narrow atomic quasi- 
momentum distribution. The Z dependence on F is very 
smooth, without the oscillations of Fig. l6la). The Niu- 



Landau-Zener prediction (dotted lines for 7). Upper two pan- 
els: Vq — 2; lower two panels: Vb = 4. 
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FIG. 5: (color online) P(t): experimental results (squares) 
and numerical solution of the Schrodinger equation describing 
the atomic cloud within the accelerated optical lattice (blue 
dashed line). The (red) continuous lines are exponential fits 
to the experimental data based on Pz(t) by Eq. (30 1, whose 
crossing with the y axis yields the value of Z. In (a) Vo = 
5.8-Erec, F = 5 and in (b) V = 1E TCC , F = 0.383. Both 
cases yield Z < 1. The slope of the exponential decay gives 
the decay rate 7. 



Raizen theory [26 , 48J includes only the two lowest energy 
bands and docs not take into account tunneling phenom- 
ena such as RET or higher excited energy bands. The 
Niu-Raizen model is thus essentially a two-state model 
for Landau-Zener coupling, neglecting resonant tunnel- 
ing effects, and averaged over all quasi-momenta in the 
entire Brillouin zone. Such a model is better suited for 
large values of Vo, when the energy bands become flat. 



B. Phase control 

To further verify that the phase <p is, indeed, the im- 
portant quantity determining the temporal evolution of 
the atomic wave function, it could be interesting to per- 
form a LZ experiment for which the atomic acceleration is 
stopped after each Bloch period for a time t halt = ir/AE, 
with AE the energy difference between the two bands, 
in order to reverse the phase of the wave function's evo- 
lution. Differences in the predicted time dependence of 
P(t) with and without this phase reversal are reported 
in Fig. [2j Even if the experimental error introduced by 
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FIG. 6: (a) Scaling plot of Z - 1 vs. of Eq. §V7}, de- 
rived from RET experimental results (squares) using a nar- 
row atomic quasi-momentum distribution. The experimental 
point at (f> = 4.8 is obtained from the data of Fig.|5|a) and the 
point at = 34 from the data of Fig. [5^b) . Full lines are the 
theoretical predictions for Vo = 1,2,3,4. The RET coupling 
yields the oscillating behavior of Z vs 4>, with the oscillation 
amplitudes increasing at lower Vo for a fixed <j>. (b) Theoret- 
ical prediction for Z — 1 in an Austin-type experiment, with 
a broad atomic quasi-momentum distribution, at Vb = 3, 3.5 
and 4 (continuous, dashed and dotted line, respectively). 



the phase imprinting could be too large to derive Z pre- 
cisely in this regime, the observation of a modified decay 
rate in the presence of a phase reversal would represent 
a direct proof that <fr is responsible for the resonances in 
the decay rate. 

The survival probability obtained in an experiment 
where after each period one halts or does not halt, with 
equal probability, represents another tool for modifying 
and testing the interference in successive Landau-Zener 
processes. The change of the decay rate by this ran- 
domization is equivalent to the change that would be 
obtained via bona fide quantum measurements, as in the 
standard formulation of the Zeno effect which was exper- 
imentally oberved in [57]. It can be demonstrated that 
the same atomic evolution is obtained by performing non- 
destructive survival probability measurements after each 
Bloch period, the quantum Zeno effect being achieved 
in the limit of very frequent measurements carried out 
within a Bloch period. 
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C. Emptying the second band 

A similar interesting experimental configuration is re- 
alized by totally eliminating the second band's occupa- 
tion after each Bloch period. This could be produced 
as in the measurement protocol used in Ref. [TS], by de- 
creasing the acceleration after each tunneling event from 
the ground band down to a small value such that the 
population in the second band tunnels to the continuum 
and is not confined anymore by the optical lattice. At 
the same time the population in the lower band does not 
tunnel to the second one, and is ready to be accelerated 
once again with the original large value. In this kind of 
setup all Landau-Zener steps in the survival probability 
as a function of time would have the same height on a 
logarithmic scale, determined by s\2 only. The phase <f> 
would then be totally irrelevant for the atomic evolution. 

D. Links with quantum field theory 

Finally, from a theoretical perspective, it would be 
of great interest to explore the links with wave-function 
renormalization effects in quantum field theory. In that 
context, the quantity Z arises from an analysis of the 
propagator (enforcing probability conservation in the 
Kallen-Lehmann representation 49, 50J ) and differs from 
unity at second order in the coupling constant. Z is 
smaller than unity for stable states, but is unconstrained 
and can become > 1 for an unstable state. There have 
been a few attempts [51H57) to analyze the quantum Zeno 
effect in the decay of elementary particles, but no exper- 
iment has been performed so far. It would be interesting 
to try and mimic these effects by making use of RET in 
BECs. This would take us into the realm of quantum 
simulations. 



larger than unity. If Z < 1, the decay can be slowed down 
(quantum Zeno effect) or enhanced (anti- or inverse-Zeno 
effect), but if Z > 1, only the quantum Zeno effect is pos- 
sible gl]. 

Our analysis of the atomic evolution in terms of succes- 
sive free evolutions and tunneling processes, with inter- 
ference in the population occupations, points out that 
Landau-Zener transitions and Stiickelberg oscillations 
[53] are two facets (one could say particular cases) of 
the very complex problem of the atomic evolution within 
the periodic potential produced by the optical lattice, in 
analogy to a previous analysis by Kling et al. [58) . 

For the shallow lattice regime, we have established a 
relationship between 7, Z and = 2tt(AE(V~o)) /F . We 
have demonstrated that the Zeno regime and resonantly 
enhanced tunneling are both controlled by the same pa- 
rameter <j) in an ultra-cold atomic cloud. The resonances 
in Z can be explained by a decay following the Landau- 
Zener probability in the first Bloch period and resonantly 
enhanced decay in the following periods. In contrast, 
the Niu-Raizen description [48J applied to describe the 
non-exponential decay of cold atoms in an optical lattice 
approximates the tunneling rate from the second to the 
third band by one complete decay. In the large Vq pa- 
rameter regime the RET resonances are not important 
and do not affect the quantum Zeno effect. 

A future experiment could involve a BEC atomic 
cloud in the presence of atomic interactions jT2l [13l l39l - 
HH [SS] [SH] ■ As verified experimentally [BD] , in this case 
the tunneling probabilities are not symmetric (Sy =/= Sji) 
and the effect of the RET resonances could be enhanced 
or suppressed with attractive or repulsive interactions. 



Appendix A: Check on the interrupted atomic 
evolution 



VI. CONCLUSION AND OUTLOOK 

In the pioneering work by Raizen et al. [55] HZ] the 
focus was on the deviations from exponential decay and 
the occurrence of the quantum Zeno effect and its inverse 
[29] [30] due to repeated measurements. In the present 
article we endeavored to go further and studied Landau- 
Zener transitions [321 133] under very different physical 
conditions, both in terms of initial state and parame- 
ters. This enabled us to use these effects as a benchtest 
for the study of wave-function renormalization effects in 
quantum mechanics. We have seen that by scrutiniz- 
ing the features of the survival probability of the wave 
function that collectively describes an ultra-cold atomic 
cloud, one can consistently define Z and extract crucial 
information on its behavior. It is remarkable that Z can 
be directly measured and that its deviation from unity 
yields directly measurable consequences on the experi- 
mental observables. In addition, as the experimental pa- 
rameters are varied, Z takes values that can be smaller or 



The dynamics of interband tunneling is discussed in 



Sec. Ill and hinges on the assumption of a free phase 
evolution over the Brillouin zone, interrupted by a very 
short tunneling event at the avoided crossing, at well de- 
fined times t = Tb(tt. + 1/2) with n £ N, as in upper 
panel of Fig. [2] and in Fig. [5jb) . To check the validity 
of this assumption we use the Hamiltonian H a which de- 
scribes the time evolution in the adiabatic (energy) basis. 
H a can be obtained by expanding the state \ip(t)) of the 
system in the time-dependent energy basis 



|^)) = E a »( t )K < )) 



and applying the Schrodinger equation idt 
with the Hamiltonian of Eq. ( 12 ) to obtain 



(a n \n) + a n d t \n)) = X, o, n E n \n) . 



(Al) 



= H\ip) 



(A2) 
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Taking the inner product with (m\ and using (m|n) = 
Smn we get 




FIG. 7: (color online) Adiabatic coupling strength c(t) de- 
fined in Eq. ( |A4[ | and normalized to maximal coupling plot- 
ted vs. time and optical lattice depth. Comparison between 
Lorentzian ansatz (upper panel) and numerical results based 
on Eq. (Ill (lower panel). The assumption of short tunneling 
events at the avoided crossings is valid for Vb ,$ 4.5 (shallow 
lattice). 



In the ideal Landau-Zener model of equation ^ and 
Ref. [3T] this yields for c(t) a Lorentzian function of time 
in a narrow time interval centered around the Tb/2 tran- 
sition time. The Lorentzian is displayed in Fig. |7ja) for 
different values of the potential depth Vq. Figure [7jb) 
shows the numerical result for c(t) in our system. The 
model discussed in Sec. IIIII ceases to be valid when c is 
large, at the border of the Brillouin zone. A comparison 
of the two plots in Fig.[7]clarines that the approximations 
used in our analysis break down for Vq ^ 4.5. 
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